The Weakly Pushed Nature of "Pulled" Fronts with a Cutoff 

Debabrata Panja and Wim van Saarloos 

Instituut-Lorentz, Universiteit Leiden, Postbus 9506, 2300 RA Leiden, The Netherlands 

The concept of pulled fronts with a cutoff e has been introduced to model the effects of discrete 
nature of the constituent particles on the asymptotic front speed in models with continuum variables 
(Pulled fronts are the fronts which propagate into an unstable state, and have an asymptotic front 
speed equal to the linear spreading speed v* of small linear perturbations around the unstable 
state). In this paper, we demonstrate that the introduction of a cutoff actually makes such pulled 
fronts weakly pushed. For the nonlinear diffusion equation with a cutoff, we show that the longest 
relaxation times that govern the convergence to the asymptotic front speed and profile, are given 
by ~ [{m + lf - l]7rVln^e, for m = 1,2, ■ ■ ■. 
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I. INTRODUCTION 

Pulled fronts are fronts that propagate into an unsta- 
ble state, for which the propagation dynamics is essen- 
tially that they are being pulled along by the growth 
and spreading of small perturbations about the unstable 
state, into which the front propagates. Concretely, this 
means that their asymptotic speed Vas is equal to the lin- 
ear spreading speed v* ofperturbations around the un- 
stable state, Vas = V* [^|-|6|. Fronts that propagate into 
an unstable state but for which Vas > v* are often termed 
"pushed" . The name stems from the intuitive idea [0J|] 
that in this regime, the dynamics in the nonlinear front 
region or the bulk region behind the front actually drives 
the front propagation: effectively it pushes the front from 
behind, and the front moves with a speed that is higher 
than the natural speed with which small perturbations 
about the unstable state spread by themselves ahead of 
the front. 

It is clear from the definition that the concept of a 
pulled front essentially pertains to a continuum formu- 
lation of the relevant dynamical variables. The linear 
spreading speed v* is defined and calculated in practice 
by considering perturbations of arbitrarily small ampli- 
tude about the unstable state of the dynamical equations; 
the value of v* then follows from an asymptotic analysis 
of the linearized dynamical equations Q. However, in all 
cases, in which one cannot ignore the fact that matter 
is made of discrete particles, one cannot perturb the un- 
stable state by any arbitrary small amount, because this 
amount must be at least one "quantum" of particle large. 

To model this discrete nature of the constituent par- 
ticles by means of a continuum equation, Brunet and 
Derrida lq| studied the nonlinear diffusion equation 



d<l) _ d^<j> 
dt dx^ 
with a cutoff e in the growth term f{(j))^ 

f{4>) = e(^-e)[0 - ^"], 

n > 1 , e.g., n = 2 or 3 . 
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Without the cutoff, i.e., for e = 0, this equation is the 
well-known nonlinear diffusion equation, which has been 
used since long as the simplest model to study front prop- 
agation into an unstable state |s|-[Tl|] . Brunet and Derrida 
B found that the asymptotic front speed Vas goes as 



(3) 
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where v* — 2 is the asymptotic speed of the correspond- 
ing pulled front of Eq. for e = 0. The above formula 
shows that the front speed v^ converges very slowly to the 
asymptotic speed v*; this illustrates that unlike pushed 
fronts, pulled fronts are very sensitive to small changes 
in the dynamics of the phase into which they propagate. 

In comparing with stochastic models of particles on a 
lattice, Brunet and Derrida associated the cutoff e with 
1 /N , where N is the average number of particles in a cor- 
relation region in the saturation phase behind the front 
Although the validity of this identification has been 
the matter of some debate, it appears that Eq. (||) with 
e = does give the proper asymptotic correction to 
the front speed even for very large N . We refer to the 
literature []6|,p^-p^ for a further discussion of the appli- 
cability of these ideas to stochastic models. 

It is intuitively clear that as soon as we introduce this 
cutoff, fronts that are pulled for e = 0, must actually 
become weakly pushed as soon as e > 0. After all, any 
perturbation around the value = does not start to 
grow until the local value crosses e, so strictly speak- 
ing, the linear spreading velocity w*(e) of arbitrarily small 
linear perturbations about the state (p — Q vanishes. As 
We > v*{e) = 0, one clearly must have a weakly pushed 
front. With this idea in mind, it is natural to address the 
convergence of the front speed to the asymptotic value, 
since it is well known that the speed of pulled fronts 
relaxes algebraically slowly to the asymptotic value v* 
[^-^, while pushed fronts normally have exponential re- 
laxation to their asymptotic speed. 
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These observations motivate us to investigate here the 
slowest relaxation modes of the stability spectrum of 
fronts for the nonlinear diffusion equation (|l|), with a 
cutoff e in the growth term (H). We calculate these re- 
laxation modes explicitly for small e, and find that the 
slowest relaxation times t„j are given by 



[(m + 1)2 - l] 



1, 2. 



(4) 



Hence, the relaxation times of the front velocity and pro- 
file approach zero as e ^ 0, but only logarithmically 
slowly. Just like the corrections to the front speed for 
practical values of e are often significant, so is the expo- 
nential relaxation, for example, for e = 10^'^, the longest 
relaxation time ti is about 4.48. Thus, while in the ab- 
sence of a cutoff the front speed is approached very slowly, 
only as 3/2t where t is the time |^-^, with a realistic 
value of e, the front speed converges relativly quickly to 
the asymptotic value. 



II. STABILITY ANALYSIS OF THE 
ASYMPTOTIC FRONT SOLUTION 

A. The Stability Operator 

The asymptotic shape of the front is a uniformly trans- 
lating front solution (j)^(x,t) which is a function of only 
the comoving coordinate = x — v^t, and which is ob- 
tained by solving the ordinary differential equation 



de 



(5) 



In carrying out the linear stability analysis of this front 
solution it is convenient to follow the standard route 
of transforming the linear eigenvalue problem into a 
Schrodinger eigenvalue problem (|J^. We consider a 
function (j){x,t), which is infinitesimally different from 
(f>e{^) = (jjeix — Vet) in the comoving frame, i.e., (/)(a;, t) — 
(j)e{x—Vet)+rj{£^, t). Upon linearizing the dynamical equa- 
tion in the comoving frame, one finds that the function 
r](x,t) = ri(^,t) obeys the following equation: 



dij drj d^f] Sf{(f)) 



(6) 



Since this equation is linear in 77, the question of sta- 
bility can be answered by studying the spectrum of the 
temporal eigenvalues. To this end, we express 7/(^, t) as 



- Et 



(7) 



which converts Eq. (g) to the following one-dimensional 
Schrodinger equation for a particle in a potential (with 
SV2m = 1): 



df2 



4 



Sf{4>) 



MO = EMO- (8) 



In Eq. d), the quantity V{0 = 



plays the role of the potential. If we now denote by 
the coordinate of the point where (jieiO = then for the 
nonlinearity (H) the potential V{S,) is easily seen to have 
the form 



(9) 



The (5-function in Eq. appears from the functional 
derivative in Eq. (0), since there is a discontinuity of 
magnitude e in at cj) — e. This discontinuity con- 
tributes an amount equal to 



dQ{(f>, 



\€i^o)\ 



S{^~^o) (10) 



to V{£,). If we combine this with the fact that |'/'£(Co)| = 
eVe, which follows immediately from the fact that one 
simply has 4>e{£,) — ee~"='^^~^°^ for ^ > one obtains 
the (5-function term in the potential given in Eq. (^. 

The form of the potential V{^) is sketched in Fig. 
|l|. Notice that is a monotonically increasing func- 

tion from e at towards the left, asymptotically reach- 
ing the value 1 as ^ — > — 00. As a result, for ^ < 
Co, V{^) also increases monotonically towards the left. 



from v;/A ~ 1 



nt^-' ~ -TrVln'e at ^ = Co-, to 
(n — vr^/ In^ e) « n as C ^ — 00. At Co, there is an attrac- 
tive 5-function potential of strength l/v^ ~ 1/2 and a 
finite step of height 1. The crucial feature for the stabil- 
ity analysis below is the fact that V{£) stays remarkably 
fiat at a value —TT^/ In^ e over a distance (Co — Ci) — U^ifl, 
and then on the left of Ci, it increases to the value « n, 
over a distance of order unity. As argued in Sec. II. B, 
this is a consequence of the nature of the solution (C) . 




width ^ I In e| " 



FIG. 1. The potential V(C) in the Schrodinger operator 
obtained in the stability analysis. 
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If there are negative eigenvalues of the above 
Schrodinger equation, then according to Eq. (0) r?(^,t) 
grows in time in the comoving frame, i.e., the front so- 
lution (peiO is unstable. On the other hand, if there are 
no negative eigenvalues, then the asymptotic front shape 
is stable, and the spectrum of the eigenvalues then de- 
termines the nature of the relaxation of (l){x, t) to the 
solution 

The full spectrum in general depends on the boundary 
conditions imposed on the eigenfunctions ipE- Here we 
consider only localized perturbations, for which we need 
to have 'ri{^, i) — + as ^ — + ±00. Because of the exponen- 
tial factor in Eq. (^, any eigenfunction ipE which van- 
ishes as ^ — > 00 is consistent with vanishing 77 towards the 
right jiot . However, for ^ — > — 00, the eigenfunctions TpE 
need to vanish exponentially fast with a sufficiently large 
exponent, so that when it is combined with the expo- 
nentially diverging term e~^'^^, they are still consistent 
with the requirement that rj vanishes for ^ — s- — cx). For 
the lowest "energy" eigenvalues, which we will investi- 
gate below, we will demonstrate that these requirements 
are obeyed. 



B. Shape of <pe{0 and the Zero Mode of the Stability 
Operator 

From the form in the potential, it is clear that the 
lowest "energy" eigenmodes, i.e., the slowest relaxation 
eigenmodes, are the ones that are confined to the bottom 
of the potential. This is the region where the nonlinear 
terms proportional to 0"^^ are negligible, and which is 
often called the "leading edge" of the front profile. For 
e <C 1, the solution of (j)e{0 this leading edge is given 

by i 

MO « ^ sin[z,e] e---? for 6 < ^ < Co - I lne| 

TT 

= ee-^'(«"«") for ^ > ■ (11) 

Here, Zi « 7r/|lne| and Zr = 1 + O(e^). The values of 
(t>e{0 "TT s-re continuous at C = ^Oi and (/'e(^o) — £• 

^ 1—1 1—1 

Although Eqs. (|l^) and (|l2|) suggest at first sight that 
the (/'e(C) has a node at ^ = 0, Eq. ( pT| ) is only valid 
in the leading edge, and 0e(^) crosses over to other be- 
havior around ^1 , which makes the front solution 0e (^) a 
mononically decreasing function of ^. The value of ^1 is 
set by the criterion that around ^1 , the nonlinear terms of 
starts to become significant, just like marks 
the point where the potential V{^) crosses over from the 
asymptotic value on the left to the bottom value. The co- 
ordinate ^1, therefore, is more or less fixed; on the other 
hand, ^0 asymptotically diverges as ~ |lne| for small e, 
making (^0 — £.1) also diverge as ~ | lne|. This is an im- 
mediate consequence of the overall exponential decay of 
4>e{0 ill £ at the leading edge. 



From the form in the potential, it is clear that the 
lowest "energy" eigenmodes, i.e., the slowest relaxation 
eigenmodes, are the ones that are confined to the bottom 
of the potential. We notice that among these modes, in- 
variably there is a zero mode of the stability operator 
that is associated with the uniformly translating front 
solution of a dynamical equation like Eq. (|l|): since 

{£,) and (j)^ + a) are solutions of Eq. for any ar- 
bitrary a, we find by expanding to first order in a that 

i^oiO — e^"^/^^^ is a solution of Eq. (||) with eigenvalue 

E = 0. From the result (Jill) for the asymptotic front so- 
lution, we then immediately get to dominant order 

tPo-sinz,^, ~^/|lne|, < C < ^o- (12) 

Furthermore, since (j)ei£,) is a monotonically decreasing 

function of ^, the solution V'o(0 = ^"""^^"jF nodeless. 

Since we know from elementary quantum mechanics that 
the nodeless eigenfunction has the lowest eigenvalue, this 
implies that all the other eigenvalues of Eq. (||) are pos- 
itive, i.e., the solution 0e(^) is stable. 

The spectrum of eigenvalues of Eq. (||) for E > 0, 
therefore, is going to determine the decay property of lo- 
calized perturbations ri{^, t) in time. We notice that for 
-E > /4 « 1, the value of the potential on the far right, 
the spectrum of eigenvalues will be continuous. However, 
we are particularly interested in the smallest eigenvalues 
E,n > for small m, since these are the eigenmodes that 
decay the slowest in time. These are the eigenvalues as- 
sociated with bound states in the potential well. 

C. Lowest Eigenmodes and Eigenvalues for e ^ 1 

As e ^ 0, the bottom well of the potential becomes 
very wide: its width diverges as |lne|. As we know 
from elementary quantum mechanics, the lowest "en- 
ergy" eigenfunctions then become essentially sine or co- 
sine waves in the potential well with small wave numbers 
k and correspondingly small "energy" eigenvalues. 

Based on the fact that the potential V{£,) on the left 
rises over length scales of order unity, we now make an 
approximation. In the limit that the bottom well is very 
wide and the k values of the bound state eigenmodes very 
small, it becomes an increasingly good and an asymptot- 
ically correct approximation to view the left wall of the 
well simply as a steep step, as sketched in Fig. H] — we 
thus approximate the potential by 

voio = n[i- em - [1 - e(e - eo)] 

In e 

-^'5(e-eo) + e(e - Co). (13) 

On the right hand side, there is an attractive delta- 
function potential at the point where the potential shows 
a step to a value close to 1. It is easy to check that the 
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prefactor of the delta- function of 1/2 is not strong enough 
to give rise to bound states with E < 0, and as a result, 
for very small values of e, the low-lying eigenmodes ap- 
proach sine waves with nodes at the position of the walls 
of the potential , 

V^m =i sin[A:„(^-^i)] . (14) 



Vo(0 



Equation (|T^) also confirms that as e ^ 0, the gap 
between the spectral lines decreases as ln~^ e, which is 
consistent with the fact that for a pulled front e = 
and the spectrum becomes gapless. Also notice that for 
the eigenvalues in Eq. dlq), the corresponding eigen- 
modes iPe{£,) decay as exp[— |^|] for ^ —oo and as 
exp[— U(:^/2] for ^ oo, which make e'"'^^'^ipEiO to go 
to zero for ^ ±oo, satisfying the boundary conditions 
discussed previously at the end of section II. A. 



width ~ I In e| 



FIG. 2. The approximate potential Vo(C) that can be used 
for calculating the low-lying modes for large widths of the 
bottom well, i.e., for small e. 

The condition that these solutions have nodes at the 
right edge of the well then yields 



(m -I- l)7r (m -|- l)7r 
Co -6 |lne| ' 



(15) 



implying that the corresponding eigenvalues are given by 



E„ 



[(to + 1)2 - l] TT^ 



0,1,2,. 



(16) 



Here, the first term between square brackets comes from 
the "kinetic energy" term fc^, while the second term 
orgininates from the value of the potential at the bot- 
tom. 

Note that for to = 0, the eigenmode sinfeg with eigen- 
value Eq is indeed the same as the zero eigenmode of Eq. 
( |l^ ) with ko = Zi, which we calculated from the shape of 
the front solution cf)^ in the leading edge. Besides veri- 
fying the consistency of our approach, this also confirms 
that there are no corrections to Eq. (jig ) for m = 0: for 
TO = it will yield an eigenvalue zero to all orders in e. 
Therefore, the smallest nonzero eigenvalue, which gov- 
erns the relaxation of the front velocity and profile to the 
asymptotic ones is Ei with relaxation time Ti given by 



rf 1 = El 



In^ e 



(17) 
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